function [ X,Y,Z ] = shift_sphere( X_o,Y_o,Z_o, centers, radius)
X = X_o*radius + centers(1);
Y = Y_o*radius + centers(2);
Z = Z_o*radius + centers(3);